Damping of dHvA oscillations and vortex-lattice disorder in the peak-effect region of strong 

type-II superconductors 
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The phenomenon of magnetic quantum oscillations in the superconducting state poses several questions 
that still defy satisfactory answers. A key controversial issue concerns the additional damping observed 
in the vortex state. Here, we show results of (iSR, dHvA, and SQUID magnetization measurements 
on borocarbide superconductors, indicating that a sharp drop observed in the dHvA amplitude just below 
H C 2 is correlated with enhanced disorder of the vortex lattice in the peak-effect region, which significantly 
enhances quasiparticle scattering by the pair potential. 

PACS numbers: 74.25.Ha, 74.25.Uv, 74.70.Dd 



The recent observations of Shubnikov-de Haas (SdH) 
and de Haas-van Alphen (dHvA) oscillations in high- 
temperature copper-oxide superconductors yj] have drawn 
increasing attention to magnetic quantum oscillations 
(MQO) as a key technique for investigating the supercon- 
ducting (SC) states of these materials at high magnetic 
fields. However, the usefulness of these studies depends 
crucially on the existence of a reliable quantitative the- 
ory for MQO in the vortex state of strong type-II super- 
conductors, a theory which does not exist today even for 
conventional type-II superconductors The demanding 
situation involved in developing such a theory may be il- 
lustrated by the extensive investigations performed on the 
nonmagnetic borocarbide superconductors YNi 2 B 2 C, and 
LuNi 2 B 2 C JIS]. Terashima et al. fl applied the field- 
modulation technique to YNi 2 B 2 C, finding a strong sup- 
pression of the dHvA amplitude just below the entrance to 
the SC state, followed by a recovery of the signal at slightly 
lower fields, and a very smooth additional damping over a 
large field range below the SC transition. The small re- 
gion of strong additional damping was found to correlate 
with that of a significant peak observed in the measured 
magnetization (the so-called peak-effect, PE). Remarkable 
persistence of the dHvA signal deep in the vortex state of 
LuNi 2 B 2 C has also been reported by Isshiki et al. 0]. 

Similarly, large additional damping of the dHvA signal, 
coinciding with the onset of the PE, was recently observed 
in YNi 2 B 2 C by employing the cantilever torque technique 
JH], with the signal persisting at significantly lower fields 
(down to 3 T). Finally, in a series of dHvA measure- 
ments, carried out very recently on LuNi 2 B 2 C by the field- 
modulation technique 0], the salient features reported for 
YNi 2 B 2 C flU have been confirmed (Fig. CO. 

The current theoretical approaches to MQO in the vor- 
tex state 0] do not provide a consistent predictive frame- 



work for a quantitative interpretation of the intriguing ex- 
perimental results. The mean-field theories based on a de- 
tailed exposition of the quasi-particle excitations obtained 
by solving the corresponding Bogoliubov-de Gennes equa- 
tions for an ordered vortex lattice 0-|9|], provide deep in- 
sight into fine features of the Landau band structure, but 
lose their transparency very quickly and become heavily 
numerical at early stages of the analysis. A simple formula 
for the additional damping ifloll . used frequently in the lit- 
erature, has been shown to be essentially valid only in the 
limiting case of a random vortex distribution 0]. The re- 
sulting additional damping rate of MQO seriously overes- 
timates the rate calculated numerically (see Fig. 8 in Ref. 
19||) for an Abrikosov vortex lattice. 

This may indicate that enhanced quasi-particle scattering 
by the SC pair-potential due to strong disorder of the vor- 
tex lattice in the PE region is responsible for the enhanced 
additional damping reported in Ref s. fllHlfJ. However, as 
proposed in Ref. yfl, an increased "phase-smearing" effect, 
i.e., broadening of Landau levels by the inhomogeneity of 
the internal magnetic field associated with displaced flux 
lines by random pinning centers, might also be responsible 
for the enhanced damping in the PE region. 

In the present paper it is shown, by means of /iSR and 
complementary SQUID magnetization measurements, that 
phase smearing is much too small to account for the en- 
hanced additional damping. It is also found that strong dis- 
order of the vortex lattice, revealed by characteristic change 
of the ^SR line shape, is correlated with the strong sup- 
pression of the dHvA signal in the PE region, whereas the 
establishment of a well-ordered vortex lattice well below 
the PE region closely follows the weak additional damp- 
ing of the dHvA signal, observed in this broad field-range. 
Both findings strongly support a scenario whereby quasi- 
particle scattering is enhanced by the disorder of the vortex 
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FIG. 1: (color online) (a) dHvA oscillation signals, for up- 
ward and downward field-sweeps respectively, measured on 
LuNi 2 B 2 C after background subtraction, and (b) the correspond- 
ing Dingle plots (triangles and circles respectively). The error 
arising from the background subtraction analysis do not exceed 
the diamater of the data points. The gray oscillatory lines in (a) 
represent the extrapolated dHvA signal, based on the normal state 
Lifshitz-Kosevich formula, whereas the solid straight line in (b) is 
the corresponding Dingle plot. The total magnetization for both 
field-sweeps (solid curves) are also shown in (b). The dash-dotted 
line in (b) represents the result of a calculation based on the ran- 
dom vortex distribution model with the zero field order parame- 
ter Ao = 4 meV and mean-field H c 2 = 8 T. Measurements were 
done at 0.5 K. 



lattice. 

High-quality single crystals of LuNi 2 B 2 C and YNi 2 B 2 C 
were prepared by the flux-growth technique at Ames Lab- 
oratory, USA flUl . and by the zone-melting method at the 
IFW Dresden, Germany fl2Tl . respectively. The SQUID 
magnetization measurements were performed at the Dres- 
den High Magnetic Field Laboratory, Germany, while the 
/iSR experiments were carried out at TRIUMF, Vancou- 
ver, Canada. The dHvA data exploited in the compara- 
tive analysis was extracted from previouse measurements 
employing the field-modulation technique on LuNi 2 B 2 C at 
T = 0.5 K |@]. A detailed discussion of their dHvA sig- 
nal in the normal state and the electronic band structure 
is given in Ref. Ill 311 . Note, that due to technical reasons 



the /iSR measurements were carried out on samples similar 
(but not identical) to the ones where the dHvA effect was 
measured. Furthermore, it is not currently feasible to carry 
out fiSR measurements below 2 K in fields above 5 T. Thus, 
in order to allow cross-correlation between dHvA and /iSR 
measurements, the magnetization was also measured using 
a SQUID magnetometer on the same samples and at the 
same temperatures as employed in the fiSR experiments. 

Transverse-field (TF) /iSR measurements up to 7 T were 
carried out on the Ml 5 muon beam line at TRIUMF us- 
ing the HiTime spectrometer, which consists of muon and 



positron detectors contained within a standard He-flow 
cryostat. The external field was directed parallel to the 
c-axis of each crystal. A fast Fourier transform (FFT) of 
the TF-juSR signal closely resembles the internal magnetic- 
field distribution P(B) UM . The measurements were typ- 
ically done by cooling the sample in a fixed field to a 
temperature between 2 and 3 K, and then measuring the 
field dependence of the TF-juSR signal. Specifically, mea- 
surements were performed on the LuNi 2 B 2 C single crys- 
tal after field-cooling in 3 and 7 T, and on YNi 2 B 2 C after 
field-cooling in 0.5 and 3 T. Several temperature-dependent 
measurements were also made. In addition, measurements 
were performed on each sample above T c at 20 K to visu- 
alize the broadening of the TF-^iSR line shape by nuclear 
moments and the field inhomogeneity of the external mag- 
net. 

A typical result from a series of systematic dHvA mea- 
surements on the flux-grown LuNi 2 B 2 C crystal |@] is 
shown in Fig. Q] for field parallel to the c-axis. The drop 
seen in the Dingle plot of the dHvA signal (originating 
from the spherical Fermi surface B13I1) just below the SC 
transition is nicely correlated with the PE seen in the mag- 
netization envelope. This feature and the weak additional 
attenuation of the dHvA signal seen below the PE region, 
essentiallyagree with the results reported for YNi 2 B 2 C 
in Refs. yt y]. One should note that the extraction of 
the dHvA signal at the PE region is somewhat ambigu- 
ous, since a reliable determination of the non-oscillating 
PE background signal is challenging. In our analysis the 
background signal was subtracted by fitting suitable poly- 
nomials to different sections of the raw-data signal. The 
Dingle data was generated using step-by-step Fourier trans- 
forms over 3 oscillations. The corresponding error bar of 
the amplitude in the PE range is about 15 percent, ensuring 
our observation of a significant additional damping there. 
For other field orientations somewhat different damping of 
the dHvA signals at the PE region was found (see e.g. Fig. 
7 in Ref. Jl). 

The upper panel of Fig. [2] shows FFTs of the measured 
/iSR signals for LuNi 2 B 2 C, after field cooling to 2.3 K in 
3 T, and the corresponding probability field distribution, 
P(B), obtained by deconvoluting each FFT curve with re- 
spect to the (practically Gaussian) reference signal. The 
onset of a large broadening of P(B) and its reversed skew- 
ness in a small field range around 6 T are apparent. The 
dramatic skewness reversal of the line shape is illustrated 
in the lower panel of Fig. [2] where the skewness parameter 
a = (AB 3 ) 1/3 1 (AB 2 } 1/2 , with AB = B - (B) |Q 
is plotted for the various distribution functions. Note, that 
a negative a is due to the presence of short-range triplet 
correlation in the absence of long-range order, characteriz- 
ing a vortex-glass phase |fl5Tl . Thus, the onset of negative 
a (around 5 T in Fig. just below the PE region (Fig. [3} 
indicates that the vortex lattice is disordered in the entire 
PE region. Remarkably, the sharp change of a is seen to 



3 



- 




7T 






6T 


y 


\ 




5.3 T 


1 




5T 


1 




4T 




3T 





-150 -100 -50 50 100 150 
B-B [G] 




FIG. 2: (color online) (upper panel) Probability field-distribution 
lines, P(B), for LuNi2B2C (solid lines) at different external mag- 
netic fields, derived by deconvoluting the FFTs (open symbols) 
of the /iSR signals. Measurements at H = 3, 4, 5, and 6 T were 
done after field-cooling to 2.3 K at 3 T, whereas at H = 5.3 and 
5.6 T were performed after degaussing at 2.7 K. The lines were 
displaced vertically for clarity and horizontally for comparison. 
The reference data used in the deconvolution were measured on 
the same sample above T c at T — 20 K. The shoulder seen on 
the 3 T FFT curve around 25 G is due to muons which missed 
the SC sample, (lower panel) Field-dependence of the skewness 
parameter, a (see text), as calculated from the field distributions 
shown in the upper panel. The line connecting the data points is 
a guide to the eyes. 



correlate with the appearance of the PE and the additional 
damping of the dHvA oscillation shown in Fig.Q] The pos- 
itive values of a near unity in a broad field range below the 
PE reflect the occurrence of a well-ordered vortex lattice in 
this region, which is seen to correlate with the weak addi- 
tional damping of the dHvA signal. One should be aware, 
however, of the kinetic nature of the PE, which is reflected 
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FIG. 3: (color online) Field dependence of the longitudinal 
magnetic moment LUN12B2C around the PE region obtained by 
SQUID magnetization measurements at T = 2.3 K (upper panel), 
and the corresponding field dependence of the /iSR field distribu- 
tion width (lower panel) calculated from the data shown in Fig. 
|2] The line connecting the data points is a guide to the eyes. 



by the history dependence of the measured /iSR data. For 
example, field-cooling in 3 T followed by an increase in 
field to 6 T results in a line shape with a negative a (see 
Fig. 2), whereas field-cooling in 7 T followed by a reduc- 
tion of the field to 6 T results in a positive a (not shown). 

Finally, our isothermal SQUID magnetization measure- 
ments, performed on both LuNi 2 B 2 C and YNi 2 B 2 C sam- 
ples, are shown in Figs.[3]and|4] respectively. The magneti- 
zation data around the PE regions are plotted together with 
the width of the TF-/iSR line shapes. The sharp maximum 
in the juSR line- width, shown for both LuNi 2 B 2 C (Fig. [5} 
and YNi 2 B 2 C (Fig.|4]l, closely follows the PE in the cor- 
responding magnetization curve. However, the magnitude 
and field range of the PE strongly depend on temperature. 
Consequently, the significant difference between TF-/xSR 
line widths observed after field-cooling in 0.5 T to 2.1 K 
and after field-cooling in 3 T to 3 K, are primarily due to 
temperature. 

Summarizing our experimental results, significant line 
broadening and reversed skewness have been observed in 
fiSR measurements, both correlated with the observation of 
a PE and an enhanced damping of the dHvA oscillations. 
The origin of this remarkable correlation is most proba- 
bly due to enhanced pinning-induced vortex-lattice disor- 
der. However, the enhanced field inhomogeneity observed 
in the PE region, approximately 30 G (maximal field dis- 
tribution width in the PE) in both LuNi 2 B 2 C and YNi 2 B 2 C 
(Figs l3l4l) . seems to be too small to attribute the additional 
damping of the dHvA amplitude to further broadening of 
the Landau levels by magnetic-field inhomogeneity. 

The following analysis establishes this conclusion: 
Imagine a charged quasiparticle (with an effective mass 
m*) moving freely in a two-dimensional spatial field pro- 
file -B(r), consisting of a large uniform part B , plus a 
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FIG. 4: (color online) (a)Field dependence of the longitudinal 
magnetic moment of YNi2B2C near the PE region obtained by 
SQUID magnetization measurements for the indicated cooldown 
conditions at 2. 1 and 3 K. (b)Field dependence of the /iSR field- 
distribution width calculated from the data obtained at 2. 1 and 3 
K. 



small nonuniform part -Bi(r) = B(r) — B . The dis- 
tribution of Bi(r) is assumed to be completely random, so 
that its ensemble average is (Bi(r)) = 0. The Landau- 
level width, tt/t r , at the Fermi energy, E F , corresponding 
to the unperturbed cyclotron frequency, uj c = eB /m*c, 
due to the inhomogeneous broadening, is given in the semi- 
classical limit, n F = E f /Huj c S> 1, by llal : = 



, where (b 2 ) = f (B 1 (r)B 1 (i f )) d 2 r' . 
For a random vortex distribution near H c2 the correlation 



length is of the order of the minimal (magnetic) length, 

a Bo = \J ch/eB Q , and (b 2 ) 
(B 2 (r)), so that: 
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The damping of the dHvA amplitude associated with the 
Landau-level broadening described above can be estimated 

from the Dingle factor: R D = exp ( — ) . This may be 



compared with the extra damping factor due to direct scat- 
tering of a quasiparticle bv the pair potential in the random 
vortex distribution limit fHOllL which takes the form: 

2 
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where A = A /Hlj c , and A is the self-consistent 
Ginzburg-Landau expression for the amplitude of the SC 
order parameter. At the PE field position: B = B PE < 
H c2 , i.e., B PE ~ 6 T, with H c2 ~ 7 T (at about 3 K), we 

2 

find A ~ 0.36n^(l — B PE /H c2 ), so that combining 
Eq. £Q with Eq. ((2| one has: 
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The Dingle plot obtained from Eq. (O, modified by ther- 
mal fluctuations |5fl, for reasonable values of the adjustable 
parameters A and H c2 J6j], is shown in Fig.[T]to agree well 
with the experimental Dingle plot in the PE region. Thus, 
Eq.© implies that the additional damping rate associated 
with the enhanced field inhomogeneity observed in the PE 
region is at least two orders of magnitude smaller than the 
observed damping rate shown in Fig. [T] 

The above estimate indicates that in order to reasonably 
account for the striking drop observed in the dHvA ampli- 
tude at the PE region one should invoke the direct influence 
of the SC pair potential on the fermionic quasiparticles un- 
der quantizing magnetic field, rather than the indirect effect 
through the magnetic-field inhomogeneity induced by the 
SC currents. Within this interpretation the observed recov- 
ery of the dHvA signal below the PE region, shown in Fig. 
[T] is explained by the reduced additional damping of the 
dHvA amplitude in the Abrikosov vortex lattice [0, 0], as 
compared to that predicted for the random vortex distribu- 
tion Jj. 
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